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Abstract 

The CMS muon system has been aligned using cosmic-ray muons collected in 2008 
and beam-halo muons from the 2008 LHC circulating beam tests. After alignment, 
the resolution of the most sensitive coordinate is 80 microns for the relative positions 
of superlayers in the same barrel chamber and 270 microns for the relative positions 
of endcap chambers in the same ring structure. The resolution on the position of 
the central barrel chambers relative to the tracker is comprised between two extreme 
estimates, 200 and 700 microns, provided by two complementary studies. With mi- 
nor modifications, the alignment procedures can be applied using muons from LHC 
collisions, leading to additional significant improvements. 



*See Appendix A for the list of collaboration members 
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1 Introduction 

The primary goal of the Compact Muon Solenoid (CMS) experiment [H is to explore parti- 
cle physics at the TeV energy scale exploiting the proton-proton collisions delivered by the 
Large Hadron Collider (LHC) H2J. The CMS experiment features a large muon tracking sys- 
tem for identifying muons and reconstructing their momenta. As with all tracking systems, 
the momentum resolution of reconstructed tracks depends in part on the alignment of detec- 
tor components in space: three translational and three rotational degrees of freedom for each 
component. The alignment procedure provides corrections to the relationships between local 
detector-bound coordinate frames and a single, global frame for all CMS tracking systems. 

The muon system consists of hundreds of independent tracking chambers mounted within the 
CMS magnetic field return yoke. Three technologies are employed: Drift Tube (DT) chambers 
on the five modular wheels of the barrel section, Cathode Strip Chambers (CSC) on the six 
endcap disks (illustrated in Figs. [T] and |2| and Resistive Plate Chambers (RPC) throughout. 
The DTs and CSCs are sufficiently precise to contribute to the momentum resolution of high- 
momentum muons (several hundred GeV/ c) assuming that these chambers are well-aligned 
relative to the CMS tracker, a one-meter radius silicon strip and pixel detector. Between the 
tracker and the muon system are electromagnetic and hadronic calorimeters (ECAL and HCAL, 
respectively) for particle identification and energy measurement, as well as the solenoid coil for 
producing an operating magnetic field strength of 3.8 T in which to measure charged-particle 
momenta (all shown in Fig. [I}. 

The CMS collaboration is developing multiple techniques to align the DT and CSC chambers 
and their internal layers. Photogrammetry and in-situ measurement devices [3] provide real- 
time monitoring of potential chamber movements on short timescales and measurements of 
degrees of freedom to which tracks are only weakly sensitive. Track-based alignment, the sub- 
ject of this paper, optimizes component positions for a given set of tracks, directly relating the 
active elements of the detectors traversed by the charged particles in a shared coordinate frame. 
Methods using tracks are employed both to align nearby components relative to one another 
and to align all muon chambers relative to the tracker. 

A challenge to track-based alignment in the CMS muon system is the presence of large quan- 
tities of material between the chambers. As a central design feature of the detector, 20-60 cm 
layers of steel are sandwiched between the chambers to concentrate the magnetic field and 
absorb beam-produced hadrons. Consequently, uncertainties in track trajectories become sig- 
nificant as muons propagate through the material, making it necessary to develop alignment 
procedures that are insensitive to scattering, even though typical deviations in the muon trajec- 
tories (3-8 mm) are large compared to the intrinsic spatial resolution (100-300 jxm). Two types 
of approaches are presented in this paper: the relative alignment of nearby structures, which 
avoids extrapolation of tracks through material but does not relate distant coordinate frames 
to each other, and the alignment using tracks reconstructed in the tracker, which allows for 
a more sophisticated treatment of propagation effects by simplifying the interdependence of 
alignment parameters. 

This paper begins with a brief overview of the geometry of the muon system and conventions 
to be used thereafter (Section^, followed by presentations of three alignment procedures: 

(a) internal alignment of layers within DT chambers using a combination of locally fitted 
track segments and survey measurements (Section |3j; 

(b) alignment of groups of overlapping CSC chambers relative to one another, using only 
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Figure 1: Quarter-view of CMS with labeled muon barrel (MB) and endcap (ME) stations. The 
steel yoke is represented by darkly shaded (red) blocks between the muon chambers. Pseudo- 
rapidities and polar angles are indicated on the top and right edges of the diagram. 
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Figure 2: Transverse slices of CMS: the central barrel wheel and representative endcap disks, 
indicating the numbering of the DT azimuthal sectors in the barrel and the CSC chamber num- 
bers in the endcaps. 
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locally fitted track segments (Section |4j; 

(c) alignment of each chamber relative to the tracker, using the tracks from the tracker, prop- 
agated to the muon system with a detailed map of the magnetic field and material distri- 
bution of CMS (Section |5j). 

Procedure (c), above, completes the alignment, relating all local coordinate frames to a shared 
frame. Its performance is greatly improved by supplying internally aligned chambers from 
procedure (a), such that only rigid-body transformations of whole chambers need to be consid- 
ered. Procedures (b) and (c) both align CSC chambers relative to one another, but in different 
ways: (b) does not need many tracks, only about 1000 per chamber, to achieve high precision, 
and (c) additionally links the chambers to the tracker. 

With the first LHC collisions, groups of CSCs will be interaligned using (b) and these rigid- 
body groups will be aligned relative to the tracker with (c). As more data become available, 
comparisons of results from (b) and (c) yield highly sensitive tests of systematic errors in (c). 

Although the ideal tracks for these procedures are muons from LHC collisions, this paper fo- 
cuses on application of the procedures using currently available data, namely cosmic rays (a 
and c) and beam-halo muons from circulating LHC beam tests in September 2008 (b). In par- 
ticular, (c) requires a magnetic field to select high-quality, high-momentum muons and concur- 
rent operation of the tracker and muon systems. The CMS Collaboration conducted a month- 
long data-taking exercise known as the Cosmic Run At Four Tesla (CRAFT) during October- 
November 2008, with the goal of commissioning the experiment for extended operation |4|. 
With all installed detector systems participating, CMS recorded 270 million cosmic-ray trig- 
gered events with the solenoid at its nominal axial field strength of 3.8 T. Due to geometric 
limitations imposed by the primarily vertical distribution of cosmic rays underground, (c) is 
performed with only a subset of DT chambers using CRAFT data, though the procedure will 
apply similarly to CSC chambers, once a large sample of inclined tracks becomes available. 

The formalism and results of each procedure are presented together. Details of the data transfer 
and the computing model which were used to implement these procedures are described in 
Ref. 0. 

2 Geometry of the Muon System and Definitions 

Muon chambers are independent, modular track detectors, each containing 6-12 measurement 
layers, sufficient to determine the position and direction of a passing muon from the intersec- 
tions of its trajectory with the layer planes ("hits"). The DT layers are oriented nearly perpen- 
dicular to lines radially projected from the beamline, and CSC layers are perpendicular to lines 
parallel with the beamline. Hits are initially expressed in a local coordinate frame (x, y, z) de- 
fined by the layers: z = is the plane of the layer and x is the more precisely measured (or the 
only measured) of the two plane coordinates. On CSC layers, the most precise measurement 
is given by cathode strips, which fan radially from the beamline [6j. Defining "local rep" as the 
curvilinear coordinate orthogonal to the strips at all points, x and local rty coincide only at the 
center of each CSC layer. 

A semi-local coordinate system for the entire chamber is defined with x, y, and z axes nominally 
parallel to the layers' axes, but with a single origin. Within this common frame, the positions 
of hits from different layers can be related to each other and combined by a linear fit into 
segments with position (x, y) and direction ( g§/ gf ) • The nominal x direction of every chamber 
is perpendicular to the beamline and radial projections from the beamline. 
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Figure 3: (a) Coordinates and residuals for a DT chamber, (b) Coordinates and alignment 
parameters for a CSC chamber. 

Residuals are differences between the predicted particle trajectories and the muon chamber 
data. Residuals can have as few as one component, from layers that measure only one dimen- 
sion in the measurement plane, and as many as four components, Ax, Ay, A4|, A-j|, from linear 



fits to all residuals in a chamber, as illustrated in Fig. 3(a) 



To compare chamber positions to each other and to the position of the tracker, consider global 
coordinates (X, Y, Z) with the origin at the symmetry center of CMS, Z axis directed anti- 
clockwise along the beamline, Y upward, and X horizontally toward the center of the LHC 
ring. Local coordinate systems are related to one another through an Alignment Integration 
Frame (AIF) that only approximates these global coordinates. It is not necessary to tightly con- 
trol the definition of the AIF, as global translations and rotations of the whole system do not 
affect any physics measurements. Specific coordinate frames, at intermediate levels between 
layers and chambers and between chambers and global, are introduced in this paper as needed. 

Corrections to positions and orientations of the layers and chambers are denoted as 5 X , 5 y , 
$z, <5cp x , S(p„, and 6^, with x, y, and z expressed in local coordinates of the layer or chamber, 
and (p x , <py, and (p z as rotations around the corresponding axis; alignment corrections are small 

enough, O(mrad), to approximately commute. For CSC chambers, A{j(p), A^|^, and S rc p are 
more appropriate than Ax, A^|, and 5 X , respectively, to take advantage of the precision of the 



cathode strips. The most sensitive CSC alignment parameters are illustrated in Fig. 3(b) 



DT layers are grouped into four-layer "superlayers," each of which measures one coordinate. 
Most DT chambers contain three superlayers, and the middle one, superlayer 2, is oriented a 
90° angle with respect to superlayers 1 and 3 to measure y positions in the chamber coordinate 
frame. DT chambers farthest from the interaction point (MB4 in Fig. [TJ contain only super- 
layers 1 and 3, and are therefore insensitive to y and CSC chambers contain six identical 
layers. 

Chambers are grouped into "stations" by their distance from the interaction region, named 
MB1-MB4 in the barrel and ME±1/1, ±1/2, ±1/3, ±2/1, ±2/2, ±3/1, ±3/2, and ±4/1 in the 
endcaps, labeled in Fig. [I] Azimuthal positions are called "sectors" in the barrel and simply 
"chamber number" in the endcaps, and these are labeled in Fig. [2] 



5 



3 Internal Alignment of DT Chambers 

The layers and superlayers of DT chambers are aligned and analyzed in two steps. First, all 
layers in the chamber are aligned using a general Millepede algorithm 0, simultaneously op- 
timizing alignment parameters and segment parameters, subject to constraints from a survey 
performed during construction. In the second step, the alignment of superlayers 1 and 3 are 
checked by fitting segments in each superlayer separately and measuring alignment errors 
from deviations between pairs of superlayer segments. 

3.1 Layer Alignment 

Measured trajectories within a single chamber do not suffer from uncertainties due to scattering 
in the steel between chambers. Segments, determined from linear fits to hits in one chamber 
only, are therefore used to determine the alignment of layers inside the chamber. The segment 
parameters depend on the layer alignment parameters within the chamber, so both must be 
resolved in a combined fit. 

Without external constraints, the combined fit does not have a unique optimum, as it is insen- 
sitive to global distortions of the chamber. Consider a shear of the chamber that translates each 
layer in x by an amount proportional to its z position and an equal shear of all segment angles: 
this new geometry has exactly the same residuals as the unsheared geometry. 

It is therefore necessary to add external survey measurements to constrain the fit; these are 
taken from two sources: (1) measurements of the wire end-pins taken during superlayer con- 
struction, and (2) photographs of reflective targets on the exterior of the superlayer ("pho- 
togrammetry"), taken during chamber construction. The averages of wire positions for all 
wires in each layer provide measurements of the X positions of the layers: RMS deviations 
from design geometry are 100 }im with 30-40 jim uncertainties. RMS deviations for the pho- 
togrammetry measurements are 200 }im in x and y, 500 jim in z, and 150 jirad in (p x , <p y , and (p z . 

The alignment of a tracking system subject to survey constraints can be expressed as the min- 
imum of an objective function with terms derived from both tracks and survey. The objective 
function is 

layers tracks . j 

X 1 = E E ( A % " A i • ^ " B < • S Pj) (^wt 2 )^ 1 (A% " Aj ■ 5i - Bj ■ 8fj 



/ 



layers targets j 

£ £ (A& - C lk • (cwy 2 ), 1 (A& - Qt • + A 



k 



2 

layers 



a) 



where 



• AXij is the residual on layer i from track j (one-dimensional and in layer coordinates 
for the layer alignment case); 

• Sj = (S x , 5 y , S z , 5tp x , 8<p y/ Sq z ) is a vector of alignment corrections for layer i; 

• Aj = ( 1 — — yy j x j^£j ~Vj ) is a 1 x 6 matrix transforming alignment 
parameter errors to residuals, dependent on the layer intersection (xj,yj) and en- 
trance angle gf . of track /; 

• 5pj is a vector of corrections to the segment parameters; 

• Bi is a matrix transforming variation of segment parameters into residuals on layer i 
(including a projection to the measured direction in the layer coordinate system); 
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3 Internal Alignment of DT Chambers 



• (chit 2 )^- 1 is the inverse of the covariance matrix of the uncertainty in the hit on layer 
i, track j (single-valued in the layer alignment case); 

• Afjt is the difference between the nominal and measured position of survey target k; 

• Cik is a matrix transforming alignment errors in layer i to corrections in the position 
of survey target k; 

• (^survey 

) k 1 is the inverse of the covariance matrix of measurement k; 



• A 



Lit 



is a Lagrange multiplier to inhibit translations and rotations of the chamber. 



Segments are modeled as straight lines because the magnetic field is negligible inside the DT 
chambers. Survey targets from wire measurements are one-dimensional (the other two coor- 
dinates are given zero inverse uncertainties in (c^vey 2 )^ 1 ) and photogrammetry targets are 
included as individual survey constraints. Wire measurements apply only to relative positions 
of layers within their parent superlayer, and photogrammetry constraints apply to all layers in 
a superlayer as a group. This is expressed in C# • Sj with terms such as 6j — [S\ + + S3 + J4] /4 
to constrain layer i in a superlayer consisting of layers 1-4. Errors in survey measurements are 
assumed to be uncorrelated. 

The x 2 is quadratic in its parameters, so it is solved by the matrix inversion method. All layers 
were aligned in 5 X , 8q x , 5<p y , and S^, using about 20 000 cosmic ray tracks per chamber. The RMS 
of the corrections is 116 }im, 58 }imd, 63 ^rad, and 49 }irad, respectively and will be studied in 
more detail in the next section. 

3.2 Test of Superlayer Alignment 

To cross-check the alignment results in superlayers 1 and 3, the track-based data are compared 
with an independent set of photogrammetry measurements. The four hits in each superlayer 
define a superlayer segment with a one-dimensional position and slope, and segments from 
different superlayers must match at a common plane. Segment residuals, Axg, are the mis- 
match of these superlayer segments, and are sensitive to relative positions of the superlayers 
as well as their internal layer alignment parameters. 

Segment residuals can be used to align superlayers 1 and 3 relative to one another without 
needing external constraints. The first term in Eq. ([TJ becomes 



tracks 

Xs 2 = E 
7 



, , dx dx \ ( s p 

V °<PyS 



^2 • ( 2 ) 



with no explicit dependence on track parameter corrections Spj. The minimization of Eq. (2| 
has a unique solution in S Xs , 5 Zs , and ^ yg . It is not susceptible to shear, for example, because the 
slopes of the segments are fixed by the layers in each superlayer. Access to additional alignment 
parameters would require knowledge of y and g| which are not available in superlayers 1 and 
3. The position of superlayer 2 cannot be determined in this way because it is the only y- 
measuring detector in a DT chamber. 



To check the partially survey-based layer alignment from Section 3.1 segment residuals before 



and after layer alignment are plotted in Fig. 4(a) Each entry in the histogram is the mean of the 



Axg distribution for a single chamber. After layer alignment, segment positions and angles are 
more consistent between superlayers 1 and 3, leading to better matching at the common plane. 

To verify the photogrammetry superlayer positions, S Xs , ^z s , and <5^ ;/S from a minimization of 
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Figure 4: (a) Distribution of means of DT segment residuals for all chambers, before alignment 
presented in a dark/ dashed histogram and after alignment in a light histogram, (b) Difference 
in z of superlayers as measured by tracks (<5 Zg ) and photogrammetry (P z ). 

Eq. Q are directly compared with the photogrammetry measurements P x , P z , and P ( p lJ of the 



relative superlayer positions. Figure 4(b) presents differences between the ^ 2 g and P z for each 
chamber. Typical values of <5 2g are 1-1.5 mm due to a glue layer not included in the design 
geometry and they agree with P z on a per-chamber basis with 580 jim accuracy. Agreement 
of Sx s and <fy yg with P x and P^, respectively have 80 fim and 50 firad accuracy. The segment 
measurements were applied to correct the internal geometry of the chambers. 



4 Alignment of CSC Chambers in Rings 

The CSC chambers in the muon endcaps were designed to overlap slightly along the edges 
of the sensitive area, such that muons passing through the narrow "overlap regions" would 
be observed by both of the neighboring chambers. This allows for a high-precision relative 
alignment of neighbors, and the relative measurements can be propagated around each ring of 
mutually overlapping chambers (illustrated in Fig. [2]). All endcap rings are internally connected 
in this way except ME±l/3. 

Although CSC chamber alignment using overlap tracks and alignment with tracks from the 
tracker (described in the next section) both determine the relative positions of CSC chambers, 
the overlap method has two advantages: (1) it achieves high precision with a small number of 
tracks, and (2) it is less prone to potential systematic errors in tracking. Since the two methods 
use different datasets in different ways, comparison between them can be used to diagnose 
systematic errors in tracks from the tracker or in their propagation to the muon endcaps. The 
disadvantage of the overlap alignment method is that it does not relate the ring coordinate 
frame to the AIF A second alignment step is therefore necessary to align the whole ring as a 
rigid body relative to tracks from the tracker. 

This section begins with a mathematical derivation of the algorithm, which is an analytic so- 
lution of Eq. ([TJ without external constraints, followed by an analysis of results from a Monte 
Carlo (MC) simulation and LHC single-beam tests. Cosmic ray results are not presented be- 
cause the procedure requires forward-pointing tracks with approximate azimuthal symmetry. 
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4 Alignment of CSC Chambers in Rings 



4.1 Ring Alignment Algorithm 

The basic strategy of CSC ring alignment is to fit segments from the same muon in each of the 
overlapping chambers independently and require them to match in position and slope for all 

chambers in the ring simultaneously. Segment parameters rep and are computed using 
only cathode strips (the high-precision rep coordinate), in a coordinate frame shared by pairs 
of chambers with z = being the plane of symmetry between them. The curvature added to 
the segments by parameterizing them with a curvilinear coordinate is negligible. Chambers 
are labeled by indices i ranging from 1 to N (N = 18 or 36, depending on the ring), with the 
convention that N + 1 refers to chamber 1. 

The general alignment objective function (see Eq. ([T]| for definitions) could be applied here, 
removing survey constraints and considering each track to consist of two "hits" with two com- 
ponents each: rcp'j and d ^ .. However, the geometry of this alignment case allows for simplifi- 
cations which make the solution analytic. Overlap tracks only connect neighboring chambers, 
i and i + 1, in such a way that residuals A(r0)j = (rf^' — 1 a ^ilil — ^ilil di "^ , 



(r(p)'j +1 and A- 



dz 



dz 



dz 



constrain A 



5 t+1 ) . Minimization of terms with the form 



A(r0)} ; 

^ dz ; 



A; 



5 l -5 



(3) 



is functionally equivalent to Eq. (JTJ with Bj ■ Spj explicitly evaluated. 

The overlap region is narrow in x and g| but wide in y, so three alignment parameters, 5 r , 
$<Py' S<pz can b e accessed from these data (illustrated in Fig. 3(b) I. The expression Aj ■ yS' — 5' +1 
expands to 

/ V 




(4) 



Furthermore, correlations between 5 



<t>y 



<5*0 1 can be replaced with an order de- 



pendence. Rotating chambers (with segments following their orientation) to make segments 
parallel (<py alignment) would change their intercepts at z = 0, but translating chambers to 
make segments continuous at z = (r<p alignment) would not change their slopes. If <py is 
aligned first and all segments are refitted, recomputing all residuals with the new geometry, 
a subsequent alignment of rep (and (p z ) does not disturb the tpy minimization. Therefore, two 
objective functions, Xi 2 an d X2 2 , can be minimized separately as long as Xi 2 is optimized first, 
and the derived geometry is used to calculate quantities in xi 2 - 

Putting all of this together, 



chambers tracks 

E E 

1 i 

chambers tracks 



A 



dz ; 



1 



S i+1 



E E AMj-f^-^j+y,^; 



A-7 1 



(5) 



where (c d ( r ^ 2 )j and (o~ r <p ) l j are one-parameter errors in the residuals. 



4.2 Monte Carlo Study 
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The sum over tracks in xi 2 can be recognized as a weighted mean and the sum over tracks in 
Xi 2 as a linear fit of A(r<p)j versus y. Evaluating them (and introducing m, a, and b as functions 
returning the weighted mean, y intercept, and slope versus y of a given dataset, respectively), 
the objective functions can be replaced with 

chambers 



1 

chambers 



(6) 



x'ia 2 = E «({A(^)} /y/ })-(^-4; 



rz'+l 



1 

chambers 



i 



where ^ differs from ^1 by a constant factor and + Xib differs from xi by a constant 



factor. To avoid rotations and twists of the whole ring, Lagrange multiplier (Yj / N) 2 is 



added to x'\r (Ei^/N) 2 is added to x^, 2 , anc ^ (Ei^L/N) 2 is added to ^fc 2 ' which favor 
solutions with minimal average corrections. Each of the three minimization problems has the 
same analytic solution, found by setting the derivatives of the objective function to zero and 
inverting the resulting N x N matrix of constants. This alignment method is applicable to any 
tracking system composed of a ring of pairwise overlapping detectors. 

The method also enables three internal cross-checks; the following must be consistent with 

zero: Yd a {{^{ r< P))> Vj})' Ei K{A(/<£)j, J/,}), and Yd m({A^|^.}). These closure tests are not 
sensitive to misalignments, but they determine whether the residuals are consistent with a 
closed loop. All three closure tests were found to be consistent with zero in these studies. 

4.2 Monte Carlo Study 

The procedure was applied to beam-halo events generated by Monte Carlo (based on a sim- 
ulation described in Ref. [8]). The simulation has approximately the same number of events 
as the 2008 LHC dataset (33 000 tracks passing through the overlap regions of CSCs), but a 
different azimuthal and radial distribution. The distribution of beam-halo events is difficult to 
predict for a new accelerator, and in fact changed frequently during the first circulating beam 
tests. Since the alignment uncertainties are statistics-limited and the population of tracks in 
each chamber is only approximately the same as in data, the simulation provides only a rough 
guide for what to expect from the data. 

From an initially misaligned detector, the procedure aligned 5* with 1.04 mrad accuracy (ini- 
tially 2 mrad), 5 r( p with 230 ]dm (initially 1000 /mi), and with 0.25 mrad (initially 1 mrad), 
determined from the RMS of differences between the aligned positions and the true positions 
of the chambers. 

The second step, aligning the internally aligned rings relative to the tracker, was studied with 
a sample of simulated muons from proton collisions. It was found that 280 ]dm ring position 
accuracy in the X-Y plane can be achieved with 10 pb _1 of collisions. This second step cannot 
be applied with beam-halo muons because they generally do not cross both the muon chambers 
and the tracker. 

4.3 Alignment Results 

In September 2008, protons circulated in the LHC, producing beam-halo muons captured by 
CMS. The majority of the beam-halo data were collected from one 9-minute fill of the anti- 
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5 Alignment of DT Chambers in a Global Coordinate System 
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Figure 5: Difference between CSC positions and their photogrammetry measurements. 
Dark/ dashed histograms are before alignment and light histograms are after alignment. 

clockwise beam. More beam-halo muons illuminated the inner rings of the negative endcap 
because they tend to follow trajectories close to the beam-line and the anti-clockwise beam tra- 
verses CMS from the negative side. All chambers in ME— 2/1 and ME— 3/1 were operational, 
so these chambers were used to test the procedure. 

Application of the alignment algorithm narrows the RMS of the A(n/>) distribution from 1.42 
to 0.98 mm. For comparison, the A(r<p) distribution of the aligned simulation has an RMS of 
1.12 mm. 

To independently verify the results, they can be compared with photogrammetry measure- 
ments. The resolution in the positions of photogrammetry targets is 300 }im [9|, which trans- 
lates into resolutions of 210 ]irs\ and 0.23 mrad for 3 rc p and 6^ respectively. Photogrammetry 
measurements cannot determine 5<p because the targets lie near the y axis of the chambers. 

Figure [5] shows the difference between chamber coordinates and photogrammetry (r<p 7G and 
$z PG )/ before and after alignment. The RMS of these distributions after alignment, which 
are convolutions of photogrammetry errors and errors in the track-based measurements, are 
340 f/m and 0.42 mrad in 5 r( p and 5^ respectively. Subtracting the photogrammetry errors in 
quadrature, one can conclude that the track-based measurement has approximately 270 f/m 
and 0.35 mrad uncertainties, in rough agreement with the prediction from simulation. In the 
absence of systematic uncertainties, several hours of similar beam-halo conditions would be 
sufficient to reduce the alignment error below the 170-200 }im intrinsic hit uncertainty for these 
chambers |6|. 



5 Alignment of DT Chambers in a Global Coordinate System 

The muon alignment procedures presented in previous sections arrange layers and sets of 
chambers in self-consistent coordinate frames, but do not relate those frames to the other sub- 
systems of CMS, in particular the tracker. In this section, a method is described to align muon 
chambers relative to the tracker by propagating the tracker tracks to the muon chambers. The 
method is equally applicable to DTs and CSCs, but cosmic rays only provide large numbers of 
tracks in the central region of the barrel. 

Muons encounter significant scattering material between every two stations. With measure- 



5.1 The Reference-Target Algorithm 
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ments expressed as chamber residuals (Ax, Ay, A^|, Ag|), this means there is a large compo- 
nent of random error in the trajectory between each measurement and the next. The residuals 
are therefore broadened beyond what would be expected from the intrinsic resolutions of the 
hits, but alignment information can be derived from the peak positions of those distributions. 

In principle, muon chamber hits could be used in the track fits to narrow the distributions 
of residuals and improve the statistical precision of the alignment results. Including those 
hits biases the tracks, coupling track parameters and alignment parameters, thereby coupling 
alignment parameters of different chambers with each other. This coupling can be resolved by 
matrix inversion [7] or by reducing the weight of the muon chamber hits and iterating ||T0|, 
but the structure of the probability distribution from scattering complicates these methods. 
Moreover, sufficient statistical precision can be achieved without muon chamber hits in the 
track fits (or equivalently muon chamber hits with negligible weight in the fit). 

This section describes muon chamber alignment using tracks determined purely by the tracker. 
Because the track parameters are not a function of the muon chamber positions, there is no mu- 
tual dependency to resolve, nor are the alignment parameters of different chambers coupled. 
However, the resulting chamber positions depend on the alignment of the tracker ||TT| : opti- 
mized muon chamber residuals does not guarantee that the combined tracker-muon system is 
globally undistorted. The first subsection below describes the algorithm with a discussion of 
propagation effects. Section pT2| presents a Monte Carlo study of the procedure, and alignment 
results are presented in Sections [53] (residuals), pT4| (cross-check) , and |5 .5 1 (momentum measure- 
ments). An alternative algorithm is under development to align muon chambers using matrix 
inversion; this is briefly described in Section [53] 




5.1 The Reference-Target Algorithm 

The "reference-target" algorithm divides the tracking volume into two regions: a "reference" 
(the tracker), in which normal track-fitting is performed, and a "target" (the muon chambers), 
in which unbiased residuals are computed from the propagated tracks. The simplicity this 
affords in the correlation matrix of alignment parameters allows more emphasis to be placed 
on the study of propagation effects. 

The shape of the residuals distribution is fitted to a parameterized Ansatz function using an 
unbinned maximum likelihood method, rather than minimizing an objective function with a 
quadratic form like Eq. jlj. This can be seen as a generalization of the standard method, because 
quadratic terms such as (Axu — Aj ■ Si) 2 / (<7hit 2 )i; are the logarithm of Gaussian likelihoods. 
Substituting a non-Gaussian Ansatz motivated by physical processes in track propagation in- 
troduces non-linearity to the derivative of the objective function which cannot be solved by 
matrix inversion. The non-linear minimization package MINUIT is used instead [12 j. 

5.1 .1 Propagation Effects 

Residuals from propagated tracks can be affected by the following effects: 

• misalignment (distributed as a ^-function in x, y, d I |, ^ because it is strictly geometric); 

• intrinsic hit resolution (negligible); 

• statistical uncertainty in the fitted track parameters (Gaussian); 

• hard Coulomb scattering off of nuclei (which has non-Gaussian tails); 

• multiple Coulomb scattering (Gaussian in the limit of many interactions); 

• background from pattern-recognition errors and noisy channels (non-Gaussian); 
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• systematic errors in magnetic field map and material budget for average energy loss 
(proportional to ql pj and q/\p\ 2 , respectively where q, pr, and \p\ are the charge, 
transverse momentum, and magnitude of the momentum of the muon); 

• systematic bias in the track source distribution (a function of the path of the track 
through the reference volume). 

All but the last two effects are included in the Ansatz. The main non-Gaussian contribution is 
from events in which the muon interacts with a small number of nuclei, such that the distribu- 
tion of deflections does not fully reach the Gaussian limit of the central limit theorem. There 
is no sharp distinction between single and multiple scattering, but it is sufficient to model the 
combined effect with a function having a Gaussian core and large tails, such as the tails of a 
Cauchy-Lorentz distribution. A smaller non-Gaussian contribution from pattern-recognition 
errors and noisy channels is also absorbed into the tails. The primary significance of the tails 
is to increase the likelihood of highly non-Gaussian residuals, and therefore stabilize the deter- 
mination of the peak. 

A magnetic field map resulting from detailed modeling and CRAFT data analysis [13] was 
used in this alignment, but the result is additionally verified by performing it separately with 
positively charged muons, negatively charged muons, and averaging the two. Magnetic field 
errors add contributions to residuals which are antisymmetric in charge, as do material budget 
errors because muons can only lose energy on average, and this deflects them in the direction 
of their curvature. Any charge-dependent effects from magnetic field map or material budget 
errors would cancel in the average. The RMS of differences with respect to not applying this 
averaging procedure are 100 ]im in S x and 0.07 mrad in S<j, , the two parameters most affected 
by magnetic field. These differences are small compared with other uncertainties. 

The possibility of systematic bias in the track source has been addressed in the context of tracker 
alignment |fl4l Il5l as weak modes in the procedure and has been studied in Ref. |TT| for the 
tracker description used here. Any unresolved global distortions in the tracker would be ex- 
tended to the muon system as well, though tracks would be guaranteed to match segments in 
the momentum range of the algorithm's application. The goal at this stage is to align muon 
chambers to the positions projected by the current best knowledge of the shape of the tracker. 



5.1 .2 The Alignment Ansatz Function 



Misalignment offsets the peak of the residuals distribution, centering (Ax, Ay, A||, A^|) at the 
values 
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(7) 



where (x,y) represents the coordinates of the track intersection with the chamber and (4|, -t|) 
the entrance angle. The above matrix is an extension of Eq. (17) in Ref. [TOl to include angular 
residuals A^f and A^|, significantly increasing sensitivity to 5$ and Sq x , respectively. 
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A Voigt distribution, or convolution of a Gaussian with a Cauchy-Lorentzian, is used to model 
the Gaussian core with large tails. The function 

r°° 1 r/2 1 /— s 2 \ 

fit; t , a, r) = j_ oo - [t _ s _ w - (r/2)2 x — exp j rfs , (8) 

has one variable t with three parameters to, cr, and T. Close to the peak, the distribution is 
approximately Gaussian (because T <C a, typically by a factor of 10), and far from the peak, the 
distribution is approximately 1/f 2 . 

The fit function for the four-dimensional residuals distribution is built from a product of four 
Voigt distributions. To account for correlations between Ax and a4|, and between Ay and A^|, 
parameters tt& x and #Ay express linear dependences between them in the fit function, and are 
allowed to float freely in the alignment fit. The correlation is simply due to the fact that an 
error in the track direction A^| introduced at a distance L upstream of the chamber causes a 
Ax m LA 4| error in position. 

Without explicitly substituting the alignment parameters, the full fit function is 



dx dy 

AXo,Al/ ,A— ,A— ,P'A*/ cr Av' Cr >dx, (T dv^Ax, r Ay , r dz ,T di// a Aj/ a At/) 
dZ0 dZQ A-j- 4x Ad7 Arf; 



/f Ax; (Ax + a Ax A^), a Ax , T Ax ^j x /^A^; A^f Q , a^, ? 

f[Ay; (Ay + a Ay A^), c7 Ay , r Ay ^) x /(^A^; A% Q , a^y, T ^ . (9) 

Substituting Axo, Ayo, A^ Q , A^| from Eq. ^ introduces the track position and entrance angle 
as four new variables and the alignment parameters as six new parameters. All parameter 
values and estimates of their uncertainties are evaluated by MINUIT, seeded by the truncated 
mean and RMS of each distribution. An example fit is shown in Fig. [6] 

The alignment fit is weighted to reduce the influence of poorly formed segments. The quality 
of a segment is quantified by x 2 1 n df/ rather than uncertainties in its parameters, so segments 
in the alignment fit are weighted by W[ = (ndf/ X 2 )i/ Hi' (ndf/ ^ 2 )/' of the segment fits. The 
objective function minimized by the alignment is x 2 = w i l°g Fj where Fj is given by Eq. Q 
for each segment i. To avoid domination of the fit by a few of the most linear segments, which 
are not necessarily from the best-determined tracks (unscattered muons), segments with the 
largest 1% of weights have been excluded. 

To resolve unmodeled non-linearities in the residuals, the procedure is applied twice, taking the 
output of the first iteration as an initial geometry for the second. No subsequent improvements 
are observed in a third iteration or beyond. 

5.1 .3 Configuration for Alignment with Cosmic-Ray Muons 

The vertical distribution of cosmic rays underground imposed geometric restrictions on the set 
of chambers that could be aligned with this algorithm. Only barrel wheels —1, 0, +1, excluding 
horizontal sectors 1 and 7 (see Fig. |2j, recorded a sufficient number of muons that also crossed 
the barrel of the tracker to perform an alignment. In addition, the fits for four nearly horizontal 
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Figure 6: Mean residuals with statistical error bars versus position (x, y) and entrance angle 
(4|, ^|) in one chamber (DT wheel 0, station 1, sector 10) from cosmic-ray data, (a) before align- 
ment and (b) after alignment. Each bin selects a narrow range of the position or entrance angle 
component under consideration, but averages over all other variables. The fit function, evalu- 
ated at zero in all other variables, is overlaid as lines before and after alignment. Asymmetries 
in the position and entrance angle distributions allow for misalignments to be manifested in 
more ways in the averaged residuals than in the projections of the best-fit function, but the fit 
properly removes these misalignments, resulting in nearly flat lines in all variables. 



5.2 Monte Carlo Study 



15 



chambers (in wheel, station, sector (-1, 2, 8), (+1, 3, 8), (-1, 1, 12), and (+1, 2, 2)) failed to 
converge, all for reasons related to the scarcity of horizontal cosmic ray muons. 

For most chambers within the restricted set, however, cosmic-ray muons are sufficiently abun- 
dant that the measurement is not statistics-limited. It is therefore possible to apply a tight set 
of track quality requirements, to control systematic errors: 

• 100 < pr < 200 GeV/ c (nearly straight tracks; the upper limit guarantees statistical 
independence from one of the cross-checks); 

• 12 out of 12 hits in DT chambers of stations 1-3, 8 out of 8 hits in DT chambers of 
station 4; 

• at least 15 hits on the tracker track, with x 2 / ndf < 10. 

The cosmic-ray period included several on-off cycles of the magnetic field, and the full field 
(3.8 T) periods were shown to result in reproducible alignment parameters with the hardware 
system [3| and track residuals. Only datasets marked as acceptable for physics with full field 
from the CMS "run registry" |4] were used. Of the 270 million cosmic-ray triggered events 
in CRAFT, the above requirements select 100000 tracks (the reduction is primarily due the 
px selection and the requirement that cosmic-ray trajectories to pass through the tracker), yet 
statistical uncertainties are typically only 120 }im in 8 X . Systematic errors, which may amount 
to several hundred microns, dominate the alignment parameter uncertainties, and hence it is 
better to select the highest-quality muons. 

DT chambers in stations 1-3 are aligned in all six degrees of freedom, but chambers of sta- 
tion 4 are only aligned in 8 X , 8$ , and 8§ z , as these are the most sensitive alignment parameters 

without Ay and residuals (recall that station 4 chambers have no y-measuring superlayer). 
Unaligned coordinates are not allowed to float in the minimization. 

5.2 Monte Carlo Study 

To test the alignment algorithm, a large sample of cosmic rays was simulated (using CMSCGEN, 
described in Ref. [16|), tracks were reconstructed with misaligned muon chambers (2 mm Gaus- 
sian smearing in x, 4 mm in y and z, 2 mrad in (p x , (p y , and (p z ), and the algorithm was applied to 
restore the original alignment parameters using only the tracks. To focus on the accuracy of the 
algorithm itself, the tracker geometry, internal muon layer geometry, magnetic field map, and 
material distribution were modeled without errors (identical in simulation and reconstruction). 
All other detector effects were realistically modeled. 

The test was performed twice, the first time with 350 000 tracks passing all selection require- 
ments and again with 100000 tracks (a subsample). The first can be considered an infinite- 
statistics limit, as statistical uncertainties from MINUIT are typically 3-4 times smaller than 
the aligned position errors (see Tables [l] and [2). The second matches the statistical precision of 
CRAFT. Distributions of chamber position errors are presented in Fig. [7[ 

5.3 Alignment Results: Residuals Distribution 

After applying the algorithm to the CRAFT cosmic-ray dataset, the four-component distribu- 
tion of residuals is found to be centered on zero to the same degree as in the Monte Carlo 
simulation. Figure [8] presents the residuals distribution of all aligned chambers, where the 
alignment makes the distribution narrower and smoother. The size of the dataset is also large 
enough to see the non-Gaussian tails in detail, and that the simulated residuals distribution 
closely matches the real one. However, the raw residuals do not provide a sensitive probe 
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Table 1: Statistical uncertainties in s imulatio n and data: uncertainties for all chambers i = 1 to 
N are summarized by presenting yj ^ ^ cr, 2 where c,- is the statistical uncertainty in one of the 
six alignment parameters below. 



sample 


S x (mm) 


S y (mm) 


S z (mm) 


% (mrad) 


Sfy (mrad) 


% (mrad) 


350 k simulation 


0.059 


0.118 


0.248 


0.170 


0.038 


0.072 


100 k simulation 


0.106 


0.210 


0.443 


0.305 


0.069 


0.129 


100 k data 


0.117 


0.243 


0.512 


0.326 


0.074 


0.146 



Table 2: RMS of differences between aligned and true positions of chambers in Monte Carlo 
simulations (distributions in Fig.[7|). 



sample 


S x (mm) 


Sy (mm) 5 Z (mm) 


8^ (mrad) 


6<p y (mrad) 


5^ (mrad) 


350 k simulation 


0.192 


0.841 0.630 


0.417 


0.095 


0.287 


100 k simulation 


0.209 


0.889 0.836 


0.497 


0.148 


0.303 




Figure 7: Differences between reconstructed and true positions of muon chambers from an 
alignment performed with 350 000 simulated cosmic-ray tracks passing selection requirements 
(only showing chambers in wheels —1, 0, +1, all sectors except 1 and 7). 
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of the alignment accuracy, because alignment corrections are typically much smaller than the 
width of the distribution. 

For a higher sensitivity, the median of the residuals distribution of each chamber is calculated 
separately, then plotted as a distribution in Fig. [9j with the RMS of the distribution presented 
in Table [3] The median is less affected by non-Gaussian tails than the mean, and it is a different 
way of achieving this insensitivity than the Voigt fits used by the algorithm. 

5.4 Alignment Results: Cross-check 

The analysis of residuals in the previous subsection provides confidence that the alignment 
algorithm is operating as designed. This section presents a test of the aligned geometry using 
a significantly different method, namely local segment fits. 

A study of alignment parameter consistency for neighboring chambers using overlapping seg- 
ments from CSCs, as described in Section |4j would provide an ideal cross-check. However, 
most cosmic rays fall on the barrel, in which chambers do not overlap one another (with the 
exception of several chambers in station 4). The local alignment quality is therefore checked 
by comparing local segments from DT chambers in neighboring stations, propagated through 
only one layer of steel to the next station. 

For each sector in a pair of neighboring stations (MB1— >MB2, MB2— >MB3, and MB3— >MB4), 
segments are linearly propagated from the inner chamber to the outer chamber and the param- 
eters x local , gf l0Cal of the propagated segment are compared with those of the segment in the 
outer chamber, yielding two residuals, Ax local , A^| local . 

Curvature from the magnetic field is not included in this propagation, but the error is corrected 

by taking advantage of the fact that contributions to Ax local and A^| local from the magnetic 
field are charge-dependent. Segments are associated with the corresponding tracker tracks to 
identify their charge and to select high transverse momentum (pj > 50 GeV/c). Residuals 
from positively charged muons and negatively charged muons are collected separately, fitted 
to Gaussian distributions to identify the peaks, and averaged without weights. Since the mo- 
mentum spectra of positively and negatively charged cosmic-ray muons are the same at this 

momentum scale, the magnetic field contributions to average Ax local and A^| local cancel, leav- 
ing only differences from misalignments. This is the same procedure as used to test sensitivity 
to the magnetic field by the reference-target algorithm. 



Figure 10 shows the results of these averages for each sector and pair of stations, before and 
after alignment, and Table |4]presents the RMS of each of the presented distributions. The distri- 
butions are convolutions of errors in the alignment and errors in the segment-matching. These 
results therefore only set an upper limit on the systematic uncertainties of the alignment itself. 
In addition, global distortions of the combined tracker and muon system are not quantified by 
this method. This results in an upper limit of 0.7 mm in 5 X (proportional to Ax) and 0.6 mrad 
in 8,h (approximately proportional to A4|) for chambers in stations 1-3, and 1.0 mm, 0.7 mrad 
in station 4. 

The uncertainty in point resolution along the line of sight of tracks is therefore bounded be- 
tween the Monte Carlo prediction (Table |2|, which includes only known propagation and de- 
tector effects, and this diagnostic, which has its own systematic uncertainties. For stations 1-3, 
the uncertainty is at best 200 and at worst 700 }im in S x . 

Muons from proton collisions, which illuminate the endcaps, will enable a local cross-check 
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Figure 8: Residuals distributions before and after alignment using CRAFT data. The full curves 
describe the aligned Monte Carlo prediction and the light dashed curves indicate the peak of 
each aligned distribution and its Gaussian approximation. 
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Figure 9: Medians of residuals distributions by chamber (one histogram entry per chamber). 
T able 3: RMS of medians of residuals distributions by chamber (distributions in Fig j9 |. 





x (mm) 


y (mm) 


af (mrad) 


% (mrad) 


Aligned MC (100 k) 


0.159 


0.172 


0.066 


0.630 


Aligned data (100 k) 


0.190 


0.166 


0.085 


0.885 


Unaligned data (100 k) 


5.667 


2.570 


1.316 


1.605 
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Figure 10: Differences in DT chamber positions and angles between pairs of stations as mea- 
sured by locally propagated segments. Dark/ dashed histograms are before alignment; light 
histograms are after alignment. 



Table 4: RMS of pairwise position and angle differences by station (distributions in Fig.[l0]l. 





Unaligned 


Aligned 


MB1 - 


■> MB2 


Ax local (j^) 


1.82 


0.68 
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. Ay local . 

A§ (mrad) 


1.68 


0.57 


MB2 - 


MB3 


Ax local ( mm ^ 


3.20 


0.69 


MB2 - 


MB3 


, A r local . 

A If (mrad) 
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with the CSC ring method of Section|4] Since the latter has a demonstrated position uncertainty 
of 270 pm, such a comparison will considerably tighten the bound on alignment resolution 
uncertainty. 

5.5 Alignment Results: Effect on Momentum Measurement 

The motivation for the alignment effort is to correct reconstructed muon momentum distribu- 
tions, so the trajectories of cosmic-ray muons are re-fitted with the new geometry to verify that 
the resolution has improved. For sensitivity to the effect of misalignments in the muon system, 
energetic cosmic rays are selected with pj > 200 GeV/ c, a sample which is independent of the 
100 < pj < 200 GeV/c tracks used to perform the alignment. Tracks are reconstructed using 
hits from the tracker and the first muon station, a simple way to optimize high-momentum 
muon resolution by increasing the effective lever arm of sagitta measurements while minimiz- 
ing bias from radiative muon showers, which become prominent at several hundred GeV/ c. It 
also focuses on the connection between the tracker and the first muon station, a pair that was 



not tested with the segments described in Section 5.4 



The top half and bottom half of the cosmic ray trajectory are reconstructed separately, split at 
the point of closest approach to the LHC beamline. Any difference in track parameters between 



the top and bottom fits is purely instrumental: in Fig. 11 the fractional difference in curvature 

(Ak/( v / 2k) 

— ( K top — ^bottom)/ (v / 2k) where k — q/ pj) is plotted before and after alignment; 
the tracker-only reconstruction is also given, for reference. Assuming that the measurements in 
the top and bottom parts of CMS are statistically independent with equal resolutions, this plot 
represents the fractional error in the curvature of tracks, which is approximately equal to the 
fractional error in its reciprocal, pj. Some global distortions correlate misalignments in the top 
and bottom parts of the detector; this test is insensitive to such modes. The fitted resolution is 
given in Table [5] 

The majority of muons in this study have a momentum close to the 200 GeV/ c threshold be- 
cause of the steeply falling distribution of cosmic rays. At 200 GeV/ c, the momentum resolu- 
tion of the combined system is not expected to significantly exceed the resolution of the tracker 
alone (see Ref. |T7), Fig. 1.5), but the biased distribution prior to alignment has been repaired 
by alignment procedure. 

5.6 The Millepede Algorithm 

A Millepede algorithm |7| is under development to align the muon chambers relative to the 
tracker, as an alternative to the reference-target algorithm. Instead of using a general non- 
linear fitting package to minimize the objective function, Millepede linearizes the problem and 
solves it with matrix inversion. A potential application of this method is to combine in a single 



fit measurements from tracker tracks and the locally fitted segments described in Section 5.4 
thereby optimizing statistical precision. Before such a generalization is used, however, it must 
be tuned to reproduce the results of the reference-target algorithm. When configured to use 
tracker tracks only, the objective function is 

chambers tracks , . j . . 

X 2 = E E Uxj-Aj-SA (W)? ( Ax) - A r IA , (10) 



i i 



where Ax = (Ax, Ay, A^|, A^|), Aj is the matrix in Eq. (7|, and (c^esidual 2 ),^ 1 is the inverse of 
the residuals covariance matrix. 



To avoid the influence of non-Gaussian tails, noisy channels, and pattern-recognition errors 
in the determination of the peak of the residuals distribution, large residuals were excluded, 
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Figure 11: Fractional curvature difference between the top and bottom parts of CMS for muons 
with pj > 200 GeV/ c (k = q/ pj). Dark/ dashed distribution is before alignment, light is after 
alignment (using the reference-target algorithm), and the open dotted distribution is tracker- 
only. 



Table 5: RMS and Gaussian core fits of the fractional curvature distributions in Fig. 11 The fits 
include all data in the central region Ak/ ( \flx) < 2 x RMS. 





Gaussian mean (%) 


Gaussian width (%) 


RMS (%) 


Unaligned 


-2.5 ± 0.6 


8.6 ± 0.5 


12.3 


Aligned 


-0.9 ± 0.3 


4.3 ± 0.2 


6.9 


Tracker-only 


0.3 ± 0.3 


4.5 ± 0.2 


6.4 
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Table 6: Accuracy of Millepede alignment in Monte Carlo simulations (350 k tracks), varying 
the selection threshold of the excluded region. The resolution attained with a threshold of 2.5 7 
reproduces that of the reference-target algorithm presented in Table |2j 



selection threshold 


S x (mm) 


5 y (mm) 


S z (mm) 


1.0 7 


0.41 


1.05 


2.39 


I.57 


0.30 


0.89 


1.46 


2.0 7 


0.26 


0.84 


0.92 


2.5 7 


0.25 


0.84 


0.78 


3.0 7 


0.25 


0.85 


0.78 



symmetrically around the peak, for each of the four residuals components. Threshold values 
for the excluded region are derived from Cauchy-Lorentzian fits. With xq and 7 being the mean 
and half-width at half maximum, only residuals between xq — 2.5 7 and Xq + 2.5 7 enter the 
sum in Eq. (To) for all four components. The optimum threshold was chosen from the study 
presented in Table |6j which lists alignment position accuracy as a function of the threshold 
value. 

6 Summary and Discussion 

This paper has presented a variety of procedures to align different parts of the muon system 
with tracks: layers in DT chambers, CSC chambers in rings, and DT chambers in a coordi- 
nate frame shared with the tracker. The available data have been fully exploited: horizontal 
beam-halo and vertical cosmic-ray muons. Through comparisons with independent data, it 
was shown that the superlayer rep resolution is 80 }im within each DT chamber, that the rep res- 
olution of CSC chambers is 270 f/m within each endcap ring, and that the DT chamber positions 
along the tracks are known with an accuracy between 200 and 700 f/m (stations 1-3, wheels — 1, 
0, +1). 

In addition, several new techniques have been introduced. The superlayer structure of DT 
chambers permits an analysis of layer geometry in a way that uses tracks alone, and there- 
fore rigorously compares the result obtained with tracks with the results from survey measure- 
ments. The overlap of CSC rings permits an analytic solution to its alignment. Non-Gaussianity 
in the physics of track propagation through the steel yoke implies a non-linear extension to the 
general alignment method. 

Techniques which will be useful for re-aligning the muon system with early LHC data have 
been tested. The favorable distribution of muons from collisions will broaden the applica- 
bility of these methods and open new opportunities for cross-checks and diagnostics, which 
ultimately will lead to a better-understood momentum resolution for high-momentum muons 
and increased discovery reach for high-energy processes. 
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